Active feedback of a Fabry-Perot cavity to the emission of a single 
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We present a detailed study of the use of Fabry-Perot (FP) cavities for the spectroscopy of single InAs 
quantum dots (QDs). We derive optimal cavity characteristics and resolution limits, and measure photolumi- 
nescence linewidths as low as 0.9 GHz. By embedding the QDs in a planar cavity, we obtain a sufficiently 
large signal to actively feed back on the length of the FP to lock to the emission of a single QD with a 
stability below 2 % of the QD lincwidth. An integration time of approximately two seconds is found to yield 
an optimum compromise between shot noise and cavity length fluctuations. © 2010 Optical Society of America 
OCIS codes: 230.5590, 120.2230. 
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1. Introduction 

A quantum dot (QD) is a structure which confines elec- 
trons in all three dimensions and behaves, as a result, like 
an artificial atom with discrete energy levels [1,2]. Self- 
assembled InAs/GaAs QDs, for example, have optical 
transition wavelengths in the near infrared, fluorescence 
linewidths of the order of 1 GHz, and observable fine 
structure. Their energy levels can be tuned by a variety 
of externally applied perturbations, including mechani- 
cal stress [3], electric fields [4,5], magnetic fields [6,7] 
and optical fields [8,9]. The workhorse for spectroscopic 
studies of QD fluorescence has long been a grating spec- 
trometer employing a cooled CCD (charge-coupled de- 
vice) detector. These devices provide high sensitivity and 
large throughput, but the resolution is in practice lim- 
ited to about 7 GHz. Thus, they are not able to resolve 
the radiative linewidth of single QDs. For this reason, re- 
searchers have recently begun to employ scanning Fabry- 
Perot (FP) cavities [5, 9] as higher-resolution probes of 
QD emission. 

FP cavities are traditionally employed in two distinct 
manners for the spectroscopy of incident radiation. In 
one case, the cavity length and thus resonance frequency 
is swept. Hence, one can sweep a narrow transmission 
passband over the frequency range of interest. In the sec- 
ond case, the cavity resonance frequency is offset from 
the optical source by approximately a half-linewidth, 
so that variations in optical frequency are converted 
into variations in transmitted power. In both cases, the 
narrowband nature of the FP transmission limits the 
throughput relative to that of a grating spectrometer; 
that is the price to be paid for higher resolution. This 
is particularly important for the case of single QD spec- 
troscopy, in which the available signal is very weak. 

In this paper, we begin by deriving the FP cavity 
linewidth that gives the optimum shot-noise limited res- 
olution for a given QD emission spectrum. Then we work 
out the resulting shot-noise limits to the uncertainty with 



which measurements of QD emission frequency can be 
made with a fixed cavity. Next, we illustrate our setup 
and show representative QD spectra, with widths of ap- 
proximately 1 GHz. Finally, we use active feedback to 
lock the length of a FP cavity to the emission of a single 
QD. The feedback loop is implemented in discrete time 
steps using a digital controller, and and we have calcu- 
lated its characteristics in terms of the integration time 
and gain parameters. The fluctuations reflect a compro- 
mise between shot noise, which favors longer time steps, 
and cavity drift, which favors smaller steps. While the 
main focus of this work is to treat the QD as a frequency 
standard to which the FP cavity is locked, the same ideas 
should apply to the inverse situation, in which external 
perturbations are used to lock the emission of a QD to 
a stable FP cavity. 
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Fig. 1. (Color online) Schematic of the experiment 
measurement setup. 
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Fig. 2. (Color online) (a) Photon count rate vs FP piezo 
voltage of a QD fluorescence. The difference between the 
two sets of (fine-structure split) peaks gives a measure 
of the cavity free spectral range, FSR ss 40 GHz. The 
smaller set of split peaks corresponds to the excitation 
of a higher order mode in the FP cavity. Inset: QD 
spectrum obtained from measuring the QD fluorescence 
in the spectrometer, without the FP cavity. The tallest 
peak in this spectrum is selected by the spectrometer 
and resolved with the FP cavity to give the data shown 
in (a), (b) High resolution spectroscopy of a single QD 
fluorescence line (dark blue points). The Lorentzian fit 
(solid green line) gives a full-width at half maximum 
(FWHM) of 1.2 GHz. For comparison, the cavity line 
is plotted on the same graph (red curve). The cavity 
has a FWHM of 0.25 GHz, implying a QD FWHM of 
0.95 GHz. 



2. Fabry-Perot analysis of QD spectra 

A FP cavity consists of two mirrors separated by a 
length L, and has a periodic transmission spectrum 
with peaks separated by the free spectral range FSR = 
c/(2L), where c is the speed of light. In the vicinity of 
a cavity resonance v cav , the peaks are approximately 
Lorentzian, 
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where T is the resonant power transmission (T < 1). 
The full-width at half maximum (FWHM) is given by 
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where the cavity finesse F is related to the fractional 
round trip power loss P r t by F = 2it / P rt . The resonant 
transmission T depends on the cavity losses and the de- 
gree to which the cavity is impedance-matched [10] and 
mode-matched [11]. 

We now consider the transmission of QD fluorescence 
with a Lorentzian spectral distribution centered at fre- 
quency v q d and having FWHM Av q d'- 



QD 



Av qd 



1 



2tt 



(^-M 2 + (^) 2 



(3) 



This expression is normalized such that the integral over 
all frequencies gives the total photon emission rate Tqd 
sent to the FP cavity. Considering a particular cavity 
resonance v cav near v q d, the cavity transmission is the 
convolution of the QD fluorescence with the FP response, 
yielding a Lorentzian in v cav centered at frequency v q d 
(see Fig. 2b): 



T(v cav ) = r 
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with FWHM, Aisdet, given by the sum of the individual 
lincwidths 

Av de t = Av qd + Av cavi (5) 

and where 

Av cav 
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is the peak detection rate. Equations (5) and (6) indi- 
cate a tradeoff between throughput (improved by large 
Av cav ) and resolution (improved by small Av cav ). 

Next, we address the situation in which the FP cavity 
length is constant and v q d varies in response to some 
externally applied perturbation. If a FP transmission 
resonance is fixed at an offset of approximately a half- 
linewidth from the center of a QD fluorescence peak, 
Wcav ~ v qd \ ps Avdet/2 (as illustrated in Fig. 2b), the 
throughput is, to lowest order, linearly related to the fre- 
quency offset Vcav — v q d- In the vicinity of this point we 
now calculate the shot- noise limited uncertainty, 5v qd , in 
a measurement of the QD offset frequency. Changes in 
v q d arc related to changes in the detected transmission 
ST by 

Shot noise imposes an uncertainty [12] 
ST = vT77, 



<r, (8) 

on a measurement of T made during an interval r. The 
corresponding limit to the measurement sensitivity is 
found by substituting equation (8) into equation (7) and 
minimizing the resulting expression for Sv q d- One finds 
that Sv q d is minimized for detunings v cav — v q d such that 
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At this bias point the measurement resolution imposed 
by shot noise is 



5v q d = 



3^3 {Av qd + Av cav f' 2 



QD ' 



(10) 



Differentiating with respect to Av cav , one finds that the 
optimal shot-noise limited resolution is obtained for 
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yielding a resolution of 
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It has implicitly been assumed in the preceding dis- 
cussion that the background (dark) count rate, T^, is 
negligible. If this is not the case, a near-optimal bias 
point is found by neglecting the shot-noise contribution 
to the noise and maximizing dT /dv q d- This bias point 
corresponds to the situation where 



r(v opt ) = -r ma * + r bk . 
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Proceeding as above, the optimal resolution is obtained 
for 

AVcav = Av qd , (14) 



with a shot-noise limited resolution of 
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3. Experimental implementation 



The FP cavity used in this work employs concave (ra- 
dius of curvature 25 mm) dielectric mirrors with an 
average reflectivity of 98 % over the wavelength inter- 
val 875 nm to 1000 nm, yielding a typical finesse of 
F w 160 near the center of the wavelength range. It 
has a length of 3.8 mm, corresponding to a free spec- 
tral range FSR = 40 GHz and a lincwidth Av cav of 
250 MHz. It was constructed by gluing both mirrors di- 
rectly to an annular piezoelectric (PZT) actuator. The 
cavity is approximately impedance-matched by virtue of 
the fact that the mirrors have equal reflectivities. Mode- 
matching was accomplished by means of a single achro- 
matic lens between the single-mode fiber used to collect 
the QD emission and the FP cavity. With this setup we 
obtained up to 60 % transmission on the resonance peak 
of the TEMqo mode of the FP. 

Our sample consists of InAs QDs embedded in a pla- 
nar distributed Bragg reflector (DBR) cavity [13]. The 
DBR cavity [14] is resonant at 920 nm and consists of 
a wavelength sized GaAs spacer at 920 nm, sandwiched 
between 15 AlAs/GaAs quarter- wavelength pairs on the 
bottom and 10 pairs on the top. The directional emission 
of the DBR cavity increases our collection efficiency, en- 
abling us to resolve QD fluorescence peaks on timescales 



as low as 1 s. We excite the QDs by pumping with a 
helium-neon laser above the GaAs band gap, resulting 
in QD emission around 920 nm. The QD density is on 
the order of 1 /j,m~ 2 and hence we can isolate single QDs 
by using a tightly focused laser beam. 

A diagram of the complete experimental apparatus 
is given in Fig. 1. The sample was maintained at ap- 
proximately 4 K in a homemade closed-cycle cryostat. 
Excitation light at 633 nm was coupled into the cryo- 
stat via fiber and focused onto the sample by means of 
a homemade cryogenic microscope objective (numerical 
aperture ~ 0.4). QD fluorescence was captured by the 
same objective and coupled back into the same fiber. 
The microscope objective was scanned over the surface 
of the sample using cryogenic positioning stages until a 
strongly emitting QD was located. Initially, the FP cav- 
ity shown in Fig. 1 was not used, and the emission spec- 
trum was measured using a grating spectrometer with 
a resolution of approximately 8 GHz. The spectrome- 
ter was used to identify single QD emission lines (inset, 
Fig. 2(a)). Next, we coupled the QD emission into the 
scanning FP cavity and then injected the transmitted 
light into the grating spectrometer in order to suppress 
light from all modes of the FP cavity but one; these 
modes could potentially transmit light from other QDs 
located within the laser excitation region. The output of 
the spectrometer was directed to a single photon count- 
ing module (SPCM). When scanning the FP cavity, a 
field-programmable gate array (FPGA) card was used 
to correlate photon arrival times with PZT voltage, and 
histograms of the number of counts vs FP cavity length 
were acquired. 

Typical QD spectra are shown in Fig. 2. The QD spec- 
trum in Fig. 2a has a fine structure splitting [15], con- 
tains 10 5 points and took about 5 minutes to acquire, 
whereas the QD spectrum in Fig. 2b is a single peak 
with 10 points and took about 1.5 minutes to acquire. 
A Lorentzian fit to Fig. 2b yields a peak width (FWHM) 
of 1.2 GHz, implying a QD fluorescence spectral width 
of 950 MHz after deconvolving the broadening from the 
FP cavity (equation (5)). The integration time was suffi- 
ciently short that we do not believe that drifts in the FP 
cavity length significantly contributed to instrumental 
broadening. Independent measurements show the radia- 
tive lifetime of similar QDs to be longer than 0.8 ns, 
corresponding to a radiative linewidth below 200 MHz. 
Hence, the emission lines appear to be broadened with 
respect to the radiative lincwidth. This phenomenon is 
usually attributed to spectral diffusion and is frequently 
observed [16-18]. 

4. Active stabilization of cavity to QD emission 

We next discuss using the signal transmitted by the FP 
cavity to lock the cavity resonance to that of the QD. 
This technique could be used in an experiment involv- 
ing a long data collection run, during which the cavity 
may otherwise drift significantly. The data acquisition 
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Fig. 3. (Color online) (a) Measured standard deviation 
of transmitted photon count rate T as a function of inte- 
gration time t, for P = 0.05, 1=1 (circles) and P = 0.1, 
1 = 1 (squares). Solid curves: Corresponding theoreti- 
cal predictions assuming power spectral density of cavity 
length fluctuations varies as uj~ 2 - 9 . Dashed curve: The- 
oretical prediction for P = 0.05, 1=1 assuming power 
spectral density of cavity length fluctuations varies as 
a;" 12 , (b) Solid points: Measured standard deviation of 
F as a function of P, for r = 0.11 s (triangles), r = 0.89 s 
(circles), and r = 3.56 s (squares); integral gain I = 1 in 
all three cases. 



can be intermittently halted and the cavity locked to 
the QD emission temporarily in order to reestablish the 
PZT offset voltage corresponding to the QD spectral fea- 
ture. Alternatively, the method may be used to lock the 
FP cavity to the side of a QD fluorescence peak so that 
the cavity functions frequency discriminator. The 
transmission can then be used to infer the response of 
the QD emission frequency to an external control param- 
eter, such as stress. Although in many applications the 
feedback would be intermittent, in the interest of sim- 
plicity, the following discussion assumes the feedback to 
be on continuously. 

In our implementation, the FP cavity is locked to the 
side of the QD emission peak. Photons are detected at 
the output of the FP cavity and counted for a integration 
time t to infer a rate T. The difference between the meas- 
ured count rate and a target count rate T° is taken as 



an error signal, and sent to a proportional-integral (PI) 
controller [19] implemented in software in the FPGA 
card. One of the FPGA analog outputs sends out a volt- 
age to control the FP cavity length in discrete time steps. 
Two sources of fluctuations are considered: Shot noise, 
with variance T/t, is present with every measurement 
of r, and the cavity is subject to environmental distur- 
bances in length 5L that change its resonance frequency 
on a slow time scale. The dynamics are controlled by 
the integration time r as well as the gain parameters P 
and /. Intuitively, a smaller integration time results in 
a more agile loop capable of responding more quickly to 
changes in the cavity length, at the expense of increased 
shot noise. 

Any feedback loop that forces a lock to the side of 
a resonance feature has the undcsircd effect of convert- 
ing amplitude noise to frequency noise. In the present 
case, amplitude variations could come from pump laser 
power fluctuations or gradual misalignment of the collec- 
tion optics and QD. In an independent experiment, we 
monitored the QD photoluminescence over time without 
the FP cavity. We found that for long integration times, 
where the cavity drift noise dominates the shot noise, 
the frequency noise that one would infer from the ob- 
served QD amplitude fluctuations was small compared 
to the cavity drifts. If laser power fluctuations are sig- 
nificant, one could split off some of the reflected pump 
power from the QD signal for normalization. Alterna- 
tively, one could render the feedback first-order insen- 
sitive to changes in signal amplitude by dithering the 
cavity, employing phase-sensitive detection, and locking 
the transmission to the center of the QD resonance [20] . 
While such an implementation would be straightforward 
with our software-based controller, we limit the dis- 
cussion to the side-lock case for simplicity. 

The dynamics of the feedback loop are calculated in 
the Appendix, where the conditions for stability are de- 
rived, as well as the fluctuations expected in response 
to measurement shot noise and cavity drifts. Assuming 
variations in the cavity length arc described by a power 
spectral density (PSD) 



Ssl(u) 



(16) 



where ui is the angular frequency and a and b are con- 
stants, the variance in the steady-state photon count rate 
is found to be 



p2 
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Here /3 is the response of the cavity transmission to 
changes in its length, P and I are the proportional and 
integral gains, and J(P,I) and Q(b,P,I) are integral 
functions whose explicit forms are given in the Appendix. 
The first term decreases with integration time r, and re- 
flects the shot noise. The second term increases with r, 
and reflects cavity drifts. 

We measured the stability of the cavity used in the 
measurements described earlier by monitoring the trans- 
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mission of a frequency-stabilized laser. The FP cav- 
ity was initially offset by about a half-linewidth from 
the laser and the transmission was monitored to infer 
changes in the cavity length. Subsequent spectral analy- 
sis showed that the data were well described by a power 
spectral density 

S S l(w) oc u- L2 nm 2 /Hz (18) 

for frequencies above 3 mHz. 

Next, we implemented the feedback loop with the same 
FP cavity, using a QD with a peak count rate of T max = 
3.5 kHz and a deconvolved linewidth Av q d ~ 1.3 GHz, 
forcing a nominal count rate of T = 1.8 kHz. The stan- 
dard deviation of the steady-state fluctuations are shown 
in Fig. 3a as a function of integration time r for two dif- 
ferent sets of gain parameters P and I. The data show the 
fluctuations to be dominated by shot noise at small in- 
tegration times, and cavity length drifts at longer times. 
Also shown (dashed line) is the theoretical result given 
by equation (17) and taking equation (18) to describe 
the cavity length fluctuations. It is apparent that the 
optimal loop performance for these gain parameters is 
obtained for integration times of approximately r = 2 s, 
whereas the theoretical prediction shows a minimum at 
approximately r = 50 s. We interpret this phenomenon 
as arising from the fact that under feedback conditions, 
the cavity is subject to additional perturbations (pri- 
marily thermal) associated with the fluctuating PZT 
drive, particularly at very low frequencies. In fact, we 
get reasonable agreement with equation (17) if we take 
Ssl(<jj) oc w~ 29 nm 2 /Hz as the PSD of the cavity while 
feedback is operational (solid curves in Fig. 3). We have 
verified this empirical PSD for a variety of sets of gain 
parameters as shown in Fig. 3b. The best case shown 
in Fig. 3a had a standard deviation 5T w 40 Hz, cor- 
responding via equation (7) to fluctuations in the offset 
of the cavity resonance from the center of the QD fluo- 
rescence peak of approximately 19 MHz, approximately 
1.5 % of the QD linewidth. An optimal choice of P, I, 
and r should allow the fluctuations to be reduced to 
about 13 MHz. Further improvement could be achieved 
by using a cavity with a somewhat larger linewidth, 
650 MHz< Av cav <1.3 GHz (equations (11) and (13)). 

An interesting alternative approach would be to con- 
struct the cavity from a material with a high thermal 
conductivity, such as copper, and actuate its length by 
means of temperature. A high level of thermal stability 
could be ensured by putting the cavity in a temperature- 
controlled environment, and the high thermal conductiv- 
ity of the copper should mitigate noise of the sort induced 
by the fluctuating PZT. 

5. Conclusion 

We have described a hybrid spectrometer for use in QD 
spectroscopy employing a FP interferometer for high 
resolution, in conjunction with a grating spectrometer 
to resolve the ambiguity associated with the periodic 



FP spectrum. For optimal noise-limited resolution, we 
show that the FP cavity should be constructed to have 
a FWHM linewidth approximately equal to the QD 
linewidth. With this hybrid spectrometer we have meas- 
ured QDs with linewidths as low as 0.9 GHz and can 
easily resolve the QD fine structure. With an appropri- 
ate cavity it should be possible to use this technique 
to resolve the frequency fluctuations contributing to the 
spectral diffusion, if they are slow enough [17]. For illus- 
trative purposes, we can apply equation (6) to the QD 
and cavity used in the previous section to find a total 
captured QD fluorescence rate of Tqd = 43.4 kHz. By 
biasing at the optimal point on a cavity with transmis- 
sion T = 1 and a linewidth Av cav = 3 GHz, one finds 
from equation (10) a shot-noise limited resolution of ap- 
proximately 10 MHz in one second of integration time, 
largely independent of Av q d as long as Av q d « Av cav . 
Assuming the QD to have a natural linewidth Av q d of 
160 MHz, corresponding to radiative decay with a life- 
time of 1 ns, the measurement resolution is equal to the 
natural linewidth in an integration time of 4 ms. 

If the FP cavity is to be used as a probe of the QD 
fluorescence frequency as a function of an external per- 
turbation, or if the FP cavity is to be used as a sensitive 
spectroscopic probe for long data collection times, it can 
be useful to use feedback to lock the cavity to the side of 
a transmission resonance. We have achieved a robust lock 
by measuring the rate at which photons are transmitted 
by the cavity and implementing a software PI controller 
to act on the cavity length so as to maintain the cav- 
ity at a fixed offset frequency. An optimal integration 
time is found to be approximately two seconds, reflect- 
ing a tradeoff between measurement shot noise and cav- 
ity length fluctuations. Finally, this work paves the way 
for the inverse experiment, in which the QD frequency is 
locked to a stabilized FP cavity using the Stark [4, 5] or 
Zeeman [6,7] effects. A particularly interesting applica- 
tion would be to lock the frequencies of two distinct QDs 
to the same FP peak, so as to create indistinguishable 
photons from distinct single-photon sources. 
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7. Appendix 

The feedback loop is described by the difference equation 
T„+i = r + Kpzt /3P' le n + lY J ^rn\+l3SL n (19) 

I m=0 J 

where we denote by e n the "error" fed to the PID con- 
troller at loop cycle n: 

e„ = T n + 5T n - r° (20) 

Here Kpzt is the response of the PZT (m/V), (3 is the 
(linearized) response of the cavity photon transmission 
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rate to changes in cavity length (Hz/m), and P' (V/Hz) 
and / are chosen in software. The quantities ST n and 
<5L„ represent the measurement shot noise and stochastic 
cavity length fluctuations, respectively. 

Subtracting T n from T„ + i, we obtain a second-order 
difference equation for T n with constant coefficients and 
a stochastic forcing function: 

r„ +1 - aT n - pr„_! = p/r° 

-P{(l + I)5T n -ST n ^} 
+0 (5L n - 5L„_x) (21) 



where we have defined 



P = -j3K PZT P' 
a = 1-P(1 + J); 



(22) 
(23) 



thus P is the (dimensionlcss) proportional gain all the 
way around the loop. Taking the z-transform [21], we 
obtain 



T(z) 



PIz 2 



(z - l)(z 2 -az-P) 
, P[l -(! + !)«] 



/3(z-l) 



<ST(*) 



(24) 
(25) 
(26) 



The condition for stability is that the poles of the char- 
acteristic equation 



<7 Z — P = 



(27) 



lie within the unit circle \z\ = 1 . After some algebra, one 
finds that the necessary conditions for stability are 



< P < 1 

< 7<2/P-2 



(28) 
(29) 



Thus both P and I are limited; too much gain will cause 
instability. 

In steady state, the discrete final-value theorem [19] 
asserts that 



lim E{TJ = lim(l- 2 - 1 )r(z) (30) 

n— >oo z—>\ 



r° 



(31) 



as desired, provided that 1^0 and P^O. The variance 
of r„ is then given by [21] 



1 

2^ 



5f? (w) du> 



(32) 



where we will express the power spectral density (PSD) 
<Sp > (w) of the discrete process (indicated by the super- 
script D) r n — r° as the sum of independent contribu- 
tions corresponding to the shot noise at each measure- 
ment and the effects of the drive terms 5T n and Sli n in 
equation (21). 



Taking the measurement shot noise to correspond to a 
stationary white noise process with zero mean and vari- 
ance equal to r°/r, the corresponding PSD is simply [21] 



s&H = r°/r. 



(33) 



The effect of the shot noise 5T on the count rate T within 
the feedback loop is governed by the transfer function 



H 4r (z) 



P[l-(1 + /)^ 

z' 2 — a z — P ' 



(34) 



with the corresponding power spectrum S^(cj) related 
to the power spectrum 5^(w) by [21] 



4 1) ( W ) = 5i 3 r ( W )|H 5r (e-)| 2 



(35) 



Similarly, the effect of the cavity length fluctuations 5L n 
on the count rate T n within the feedback loop are gov- 
erned by 

h sl(z) = -5 — 5 (36) 



az-P 



(2) * 

with the corresponding power spectrum 5p (oj) given by 



4 2) M = ^M|H, L ( e -)| 



(37) 



In fact, the cavity length varies during the integration 
time r, so that Sf L {u>) corresponds to the moving av- 
erage, sampled at times nr, of a continuous stochastic 
process with power spectrum Ssl(w)- The power spec- 
trum of the discrete process 8"L n is given by [21] 



1 



. 9 .u + 2irn.^ .uj + 2-kti. 
^ smc 2 ( )S 5L ( ) 



where the term containing the sine function arises from 
the moving average and the terms in the sum with n^O 
reflect aliasing from undersampling. 

For simplicity, we now neglect the terms arising from 
aliasing, and in addition assume a form 



Ssl(u) = a\u)\ 



(38) 



for the power spectral distribution of the (continuous) 
cavity length fluctuations. Combining the previous re- 
sults, the variance of T n can then be expressed as 



4 = - I sir M + 4 1} M + 4 2) H dw 

n Jo 

r° r p 2 ~\ nR 2 

-\l + —J(P i I)\ + ^T b - 1 Q(b,P,I) (39) 



where 



J(P,/) = 

2 + 21 + 1 2 -2(1 + 7) cos(w) 
1 + P 2 + a 2 - 2cr(l + P) cos(w) + 2Pcos(2cj) 

Q(b,P,I) = 



(40) 

diu, 

(41) 
(42) 
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2[1 - cos(cj)] sinc 2 (w/2) uj~ b 
1 + P 2 + a 2 - 2<r(l - P) cos(o;) - 2Pcos(2o;) 



(43) 

J(P, I) can be evaluated analytically. The integral 
Q(b, P, I) is finite provided b < 3. 

The first term in equation (39), reflecting the shot 
noise, can be made smaller by choosing a longer integra- 
tion time t and smaller values P and / for the feedback. 
The second term, reflecting cavity drifts, can be made 
smaller by choosing a smaller integration time r. It di- 
verges for / approaching either zero or the maximum 
allowed value 2/P — 2. 
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